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ABSTRACT 

We present Chandra and VLA observations of two galaxy clusters, Abell 160 and Abell 2462, 
whose brightest cluster galaxies (BCGs) host wide angle tailed radio galaxies (WATs). We 
search for evidence of interactions between the radio emission and the hot, X-ray emitting 
gas, and test various jet termination models. We find that both clusters have cool BCGs at the 
cluster centre, and that the scale of these cores (~30-40 kpc for both sources) is of approx- 
imately the same scale as the length of the radio jets. For both sources, the jet flaring point 
is coincident with a steepening in the host cluster's temperature gradient, and similar results 
are found for 3C465 and Hydra A. However, none of the published models of WAT formation 
offer a satisfactory explanation as to why this may be the case. Therefore it is unclear what 
causes the sudden transition between the jet and the plume. Without accurate modelling, we 
cannot ascertain whether the steepening of the temperature gradient is the main cause of the 
transition, or merely a tracer of an underlying process. 

Key words: galaxies: active - galaxies: clusters: individual Abell 160 - galaxies: clusters: 
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1 INTRODUCTION 

Wide angle tailed radio galaxies (WATs) are objects located at or 
near cluster centres, with long plumes that are often bent, and are 
generally hosted by the dominant cluster galaxy. However, not all 
radio galaxies with bent plum es mee t our definition of a WAT. We 
follow the definition o f lLeahvNl993l) . in which WATs initially have 
well collimated, kpc-scale jets which suddenly flare into diffuse 
plumes, which may be significantly bent. 

WATs have tradition ally been classified as a subset of Fanaroff 
and Riley type I galaxies l Fanar off" & Rilevl 19741) because of their 
large-scale structure. Unlike classical FRIs, which have jets with 
very wide opening angles, WATs have narrow collimated jets on 
small scales (tens of kpc) that resemble those of FRII sources. 
The disruption of WAT jets occurs at the base of the plumes, rather 
than at the end of the sourc es, as is the case for FRII sources 
lHardcastle & Sakelliouteool . 

The large scale plumes in WATs can be shaped by the inter- 
actions with their environments. For example, plumes can be bent 
backwards due to the relative motion of the host galaxy through the 
intra-cluster medium (ICM) iBegelman et alll979h . The galaxy is 
not at rest with respect to the ICM, as one would expect of the 



brightest cluster galaxy (BCG) in a v irialize d cluster, but rather has 
some (small) velocity. iGomez et alJ 1 19971) interpret this as a sign 
of some cluster-cluster interaction or merger, which disturbs the 
cluster potential. However, many WAT sources are not significantly 
bent, so that it has been argued that the plume bending is not related 
to the characteristic jet flaring. For instance, 0110+152 exhibits 
sudden jet flaring, but has relatively straight plu mes; whereas other 
sourc es, such as 3C75 (see, for example H ardcastle & Sakellioul 
2004), exhibit large scale bending, but not the sudden jet flaring. 
Attempts to explain t he large scale b ending ha ve been ex plored 
in det ail elsewhere (eg lPinknev et alll994l and lO'Donoghue et alJ 
1993). In this paper, we address the possible explanations for the 
sudden jet flaring. 

It is thought that the denser cluster environment is responsible 
for the characteristic jet-plume transition of WATs, with the clus- 
ter properties determining the location of the base of the plume; 
WATs are found exclusively in cluster environments, whilst FRIIs 
of similar radio powers avoid such environments. X-rays are the 
best source of information regarding the state of the hot ICM, so 
using long exposure, high resolution observations with Chandra , 
we can probe the state of the ICM, and determine which if any of 
the models proposed to explain the sudden jet flaring fit the ob- 
served data. 
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In this paper, we present Chandra observations of the two 
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Table 1. Cluster and WAT characteristics 





Cluster centre coordinates 












Radio core 


coordinates 


Abell Cluster 


02000 


#2000 


z 


scale 


D L 


R500 


Radio source name 


«2000 


#2000 










arcsec/ kpc 


(Mpc) 


(kpc) 








Abell 160 


01 13 03.9 a 


+15 29 43.9 


0.0447 


0.88 


180 


870 


0110+152 


01 12 59.6 


+15 29 28.7 


Abell 2462 


22 39 05. 2 6 


-17 19 53.6 


0.0737 


1.57 


320 


930 


2236-176 


22 39 11.4 


-17 20 27.3 



NOTES: 

a: coordinates taken fromlAcreman^t al| |2005 ). 
b: coordinates taken from lG6mez^t^nn99^lK 

Dl is the luminosity distance to each cluster. R500 is the radius within which the mean density of the cluster is equal to 500 times the critical 
density of the universe. 



nearby galaxy clusters Abell 160, and Abell 2462, (for details of 
cluster properties see Table0. The central cluster galaxies of both 
these clusters are WAT hosts. Abell 160 is a poor cluster (rich- 
ness class 0), and was first observed with the IPC instrument on 
board the Einstein satellite. This showed emission from the cD 
galaxy, as well as some eviden ce of extended emission from the 
cluster Ijones & Formanlll999l) . X-ray observations with ROSAT 
showed extended cluster emission that was fairly uniform in struc- 
ture iDrake et alj|2000h . However, the temperature of the cluster 
has not been previously measured by either Einstein or ROSAT . 
Abell 160 ho sts the WAT 0110+ 152. Abell 2462 is another poor 
Abell cluster. iGomez et all Il997t) working with ROSAT data sug- 
gested that the X-ray structure was somewhat irregular, and that 
some elongation existed from the centre of the cluster towards the 
east. They found a temperature for Abell 2462 of l.Slto'i keV. 
Abell 2462 hosts the WAT 2236-176. 

In Section|2| we present the new Chandra observations, and 
in Section[3] the radio observations. We describe the X-ray spatial 
and spectral analysis in Section [4] In Section [5] we discuss our 
results and the possible explanations in the light of current models 
(discussed above) for jet disruption and termination. 

We assume Ho=70 km s _1 Mpc -1 and use J2000 co- 
ordinates throughout. 



2 X-RAY OBSERVATIONS 

The clusters Abell 160 and Abell 2462 were observed with Chan- 
dra on 2002 October 18 and 2002 November 19, for ~60 and 
~40 ksec respectively. During both observations, the detectors 
were operating in the VFAINT mode. The prime instrument for 
the Abell 160 observation was ACIS-I, whilst the pointing for 
Abell 2462 was centred on the S3 chip of ACIS-S. Information on 
the Chandra observations is given in Table|2| 

The data were reprocessed using CIAO version 2.3 and CALDB 
2.23, to use the VFAINT mode cleaning techniques. Additionally, 
the data were corrected for the ACIS charge transfer inefficiency 
(CTI) using the ACIS_PROCESS_EVENTS script. 

Light curves in the energy range 0.3-10.0 keV were created 
for each data set, in regions outside the cores of both clusters, show- 
ing that neither of the observations were severely contaminated by 
background flares. Subsequently, they were subjected to a 3a clip- 
ping, using the usual CIAO script, ANALYZEXT.SL. The unfiltered 
and corrected (for periods of high background contamination) ex- 
posure times are shown in Table|2| 

We used the CIAO script MERGE_ALL to create exposure maps, 
at an energy of 1.5 keV , for both clusters. For any subsequent spa- 



Table 2. Chandra observations 



Source 


Detector 11 


Unfiltered 


Corrected 






Exposure 


Exposure 






(ks) 


(ks) 


Abell 160 


ACIS-I 


58.50 


57.49 


Abell 2462 


ACIS-S(S3) 


39.24 


38.24 



NOTES: 

a: 'Detector' refers to the chip and/or chip array that was 
used as the prime instrument for each observation. 



tial analysis (e.g. Section 4.1), we divided the cluster images by the 
exposure maps to correct for exposure effects. 

A blank-sky background file was also created for Abell 160, 
by following the relevant CIAO thread. This was used to obtain the 
radial background profiles in the subsequent analysis of Abell 160 
(Section 14.31 . as the cluster fills the ACIS-I field of view, and so 
finding areas devoid of cluster emission was difficult. The blank sky 
files were subjected to the same cleaning and reprocessing methods 
as the data events file. 

Raw count images of both clusters in the 0.3-5.0 keV energy 
range are shown in Fig.Q The pixel size is 1.96 by 1.96 arcsec. We 
note that on large scales, the emission appears to be fairly symmet- 
ric and uniform. The bright, extended source in Abell 160 at 02000= 
l h 13 m 15 s <$20oo=+15°30'57" is labelled 'Source 1' in Fig. □ and 
we investigate its properties in Section l4.2.3l In Fig.|2| we present a 
smoothed X-ray image of the cluster core, in the same energy range 
as the raw counts image. The X-ray images have been smoothed 
with a Gaussian beam of width 3 arcsec and are overlaid with the 
radio contours described in Section|3| 



3 RADIO OBSERVATIONS 

VLA data at 1 .4 GHz were obtained from the VLA archive for 
Abell 160. The data were taken in the A an d C array co nfigurations 
and were originally presented in O'Donoghue et al. (1990), to- 
gether with relevant polarization maps. The archival data were cal- 
ibrated using 3C48 and 0019-000 before CLEAN maps were made 
using the AIPS task IMAGR. The Abell 2462 radio data were orig- 
inally presented in Hardcastle & Sakellioul 120041) . together with 
polarization maps. The radio source was observed with the VLA 
at 8.4 GHz in all four array configurations, and the data were pro- 
cessed as described in Hard castle & Sakellioul <2004 . 

We present the radio data in Fig. [2] as contours overlaid on 
X-ray images of the central region of each cluster. 
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Figure 1. Chandra images in the energy range 0.3-5.0 keV of Abell 160 (left panel) and Abell 2462 (right panel). The position of the BCG in each cluster is 
marked with a cross. 
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Figure 2. X-ray images of the central regions of the clusters, with the radio contours overlaid. Abell 160 is on the left and Abell 2462 on the right. The 
resolution of the radio images is 6.4 by 6.4 arcsec for 01 10+152 and 0.70 by 0.52 arcsec for 2236-176. The lowest contours are at 0.94 and 0.18 mJybeam - 1 
for 01 10+152 and 2236-176 respectively, and contours are spaced at 1,2,4,8,16,32 and 64 times the lowest contour values. 



Both sets of data show narrow, well-collimated jets that flare 
into diffuse straight plumes. Additionally, the eastern lobe in 2236- 
176 bends suddenly through 90 degrees, 25 arcsec from the jet flar- 
ing point. Details of the radio source dimensions are given in Ta- 
ble0 



4 X-RAY DATA ANALYSIS 
4.1 Spatial analysis 

To make an initial estimate of the size of the extended galactic 
emission, and to map the transition from the galactic to the clus- 
ter atmosphere, binned, exposure corrected images (as described in 



Table 3. Radio jet properties 



Source 


Jet Direction 


Jet Length 






(arcsec) 


0110+152 


NW 


40 




SE 


40 


2236-176 


NE 


28 




SW 


28 



NOTES: 

The jet lengths are taken from 
lHardcastle & SakellioJ <2004l) 
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Section[2j were used to extract the 1-D surface brightness profiles 
for both data sets. Point sources were excluded, and for Abell 160, 
the group (see Section |5J was also excluded. Both profiles were 
centred on the radio core, and adaptively binned so that each bin 
had a minimum signal to noise ratio of 3. 

Profiles were extracted out to 5 (~260 kpc) and 3 arcmin 
(~280 kpc) (~ 0.3 Rsoo, see Table0 for Abell 160 and Abell 2462 
respectively. Backgrounds were annuli centred on the radio cores, 
from 5 to 7.7 arcmin for Abell 160 and from 3 to 4 arcmin for 
Abell 2462. The profiles were fitted with a composite model com- 
prising of three components - two (3 models for the extended galac- 
tic and cluster emission, and a point spread function (PSF) to model 
the emi ssion from the active core. The (3 models are standard (3 
models (Cavaliere & Fusco-Femiano 1976), defined by: 

-3/9+0.5 



Table 4. Best-fitting parameter values for the surface brightness fits for 
Abell 160 and Abell 2462. Errors given are la for two interesting parame- 
ters. 



E(r) 



1 + 



(1) 



where £ (r) is the surface brightness at a distance r from the cen- 
tre of the distribution, So is the central surface brightness, and 
ro is the core rad ius of the dist ribution. The PSF parametrization 
of 0331+39 from lWorrall et alJ <200ll) was used for both sources, 
since the spectrum of this source was similar to those of both 
Abell 160 and 2462 (see Section l4~2l for details of how we extracted 
and fitted the spectra). 

This composite model was fitted to the surface brightness data, 
leaving the normalizations, core radii, and (3 values of the two (3 
components, and the normalization of the PSF component free to 
vary. 

The profiles and best-fitting models are presented in Fig. [3] 
and the best-fitting parameters are listed in Table |4| We find that 
for both sources the composite model fits the data well. In both 
cases, one component is a steep /3-model fitting the galactic emis- 
sion on scales of 10-20 arcsec, while the second flatter component 
fits the extended cluster-scale emission. For clusters, (3 values of 
~0.6 (Sarazin 1988) are expected, and our data agree well with 
this. 



4.2 Spectral analysis 

To further investigate the environment in which the jets propagate, 
and how the environment may relate to the causes of the jet flaring, 
we investigated the spectral properties of the WAT host clusters 
(Section f4.2.1> and the WAT host galaxies fSection f4.2.2t . 



4.2.1 The WAT host clusters 

We accumulated counts in circular regions centred on the BCGs 
out to 5 and 3 arcmin for Abell 160 and Abell 2462 respectively. 
The background regions were concentric circular annuli extending 
from 5-7 arcmin for Abell 160 and 3-4 arcmin for Abell 2462. 
From Section 14.11 it was apparent that emission from the WAT 
host galaxy is dominant out to 10 and 5 arcsec for Abell 160 and 
Abell 2462 respectively. The emission from the BCGs was masked 
out by excluding circular regions, centred on the BGCs, of radius 
20 and 10 arcsec (Abell 160 and Abell 2462 respectively). 

Spectra were extracted using the ACISSPEC script, and the 
weighted response files were corrected for the degradation in the 
quantum efficiency of the ACIS chips, using standard C1AO tools. 
We binned the spectra such that each energy bin contained a mini- 
mum of 40 counts after background subtraction. 

The spectra were fitted in the energy range 0.5-5.0 keV with 



Source 


Component 


Parameter 




A 160 


PSF 


Norm. 

(counts/arcsec 2 ) 


4.0 ±0.01 




Galaxy 


r ga i 


Q o + i.5 
°-°-0.9 




(3 Model 


(arcsec) 






(3gal 


0.96±oi 






Norm. 


3 45+0- 7 
°-*°-0.8 






(counts/arcsec 2 ) 




Cluster 




191±g 




13 Model 


(arcsec) 






(3clus 

Norm. 

(counts/arcsec 2 ) 


n tra+0.06 
U.OO_o .05 
n 16 +0.004 
u ' 1D -0.004 




X 2 (d-o.f) 




391(313) 


A 2462 


PSF 


Norm. 

(counts/arcsec 2 ) 


20.0 ±0.01 




Galaxy 


r( Vz 


2.3 ±0.3 




(3 Model 


(arcsec) 

(3gal 

Norm. 

(counts/arcsec 2 ) 


0.74 ±0.1 

5-6l 2 : 5 




Cluster 


r c!us 


165±J 




(3 Model 


(arcsec) 

(3clus 

Norm. 

(counts/arcsec 2 ) 


67+ 01 
u - D '-o.oi 

r, i n + 0.001 

u - lo -0.001 




X 2 (d-o.f) 




90 (82) 



single MeKaL models modified by the line-of-sight hydrogen col- 
umn (TVh) for each cluster. The temperatures, metal abundances, 
(Z), and JVh were left free to fit the data resulting in the best- 
fitting values presented in Table [5] No previous X-ray tempera- 
ture measurement has been made for Abell 160, but the results ap- 
pear to b e in good agreement f or othe r X-ray observations of poor 
clusters l Sanders on & Ponmanll2003l) . The temperature found for 
Abell 262 is slightly higher than that found from previous work by 
iGomez etall fiW/). but they fitted a larger cluster area, using sig- 
nificantly fewer counts, which resulted in a poorer fit to their data. 
The A^h values that we measure are consistent with Galactic values 
(4.3 and 3.1xl0 20 cm ~ 3 for Abell 160 and Abell 2462 respec- 
tively). The values of Z are in agreement with published values of 
Z in poor clusters. The spectra are presented in Fig.[4]and Table|5| 



4.2.2 The WAT host galaxies 

We extracted spectra, as described in Section l4.2.1l from circular 
regions, centred on the WAT host galaxies, of radius 20 and 10 arc- 
sec for Abell 160 and Abell 2462 respectively. Background regions 
were concentric annuli from 20-50 arcsec for Abell 160 and 10- 
40 arcsec for Abell 2462. 

The galaxy spectra were fitted in XSPEC with a composite 
model consisting of a MeKaL model and a power law modi- 
fied by the line-of-sight hydrogen column. Energy ranges of 0.5- 
5.0 keV were used for both BCGs. The spectra and best-fitting pa- 
rameters for both BCGs are shown in Fig. [5] and Table Irrespec- 
tively. The index of the power law model, the temperature of the 
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Distance from BCG/arcsec Distance from BCG/arcsee 



Figure 3. The 1-D radial surface brightness profiles for both clusters - figure (a) is the profile for Abell 160, and figure b for Abell 2462. The solid line is the 
composite model, the dashed line indicates the contribution from the PSF (modelled as described in the text), the dot-dash line is the contribution from the 
steeper (3 model for the galaxy, and the dotted line is the contribution modelled by the second fi model. The dashed vertical line shows where the radio jets 
flare. 




channel energy (keV) channel energy (keV) 

(a) Abell 160 (b) Abell 2462 



Figure 4. The spectra of the clusters Abell 160 and Abell 2462. Both spectra were fitted with a MeKaL model. The emission from the BCG and AGN has 
been removed. The lower panel of each plot shows the residuals of the best-fitting model. 



MeKaL model, and the normalizations of both models were al- 
lowed to vary. All other parameters were fixed - absorption to the 
Galactic value (given in Table[6j, and metal abundance to O.3Z0. 

Our values for g alaxy tempe r atures appear to be similar 
to values published by Heinz et al. 1 2002 ) for the cD galaxy in 
Abell 4059 and bv iTakizawa et alJ 12003 i for the cD galaxy in 
Abell 3112, which are both radio galaxy hosts. 



4.2.3 The bright, extended source in Abell 160 - 'source 1 ' 

In Abell 160 there appears to be a bright region of emission near the 
cluster redshift (see Section^. To investigate what this could be, a 
spectrum of the region was extracted using a circular source region 



of radius 0.5 arcmin (26 kpc) and a local background in the form of 
an annulus extending from 0.5 to 1.0 arcmin (53 kpc) both centred 
on the coordinates stated in Section [2] We extracted the spectrum 
in a similar way to that described in Section l4.2.1l 

This spectrum was fitted with a MeKaL model, modified by 
absorption from the line-of-sight hydrogen column. A^h was fixed 
to the Galactic value, and the abundance to O.3Z0. Allowing the 
temperature and normalization of the MeKaL model to vary, we 
obtained a temperature of 0.52+Q'Qg keV , significantly lower than 
the temperature of the surrounding ICM. We present the spectrum 
of source 1 in Fig. [6] 

From the DSS optical data, we see that most of the X-r ay emis- 
sion is associated with a large elliptical, at 2=0.0484 IPinknev et alJ 
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Figure 5. The power law and MeKaL fits to the spectra of the BCGs in Abell 160 and Abell 2462. The solid line shows the best-fitting model, and the 
residuals to the fit are shown below the spectrum. 



Table 5. The best-fitting parameters for the clusters, fitting a MeKaL 
model, and taking into account absorption from the Galaxy. 



Source 


Component 


Parameter 




A 160 


MEKAL 


kT ( keV ) 


2.7 ± 0.2 






Z(Z ) 


0.4±0.1 






Lx 


4.0±0.02 






xlO 42 erg s - 1 








N H (XlO 20 cm " 3 ) 


4.4 ±0.1 




X 2 (d.o.f) 




177 (161) 


A 2462 


MEKAL 


kT ( keV ) 


2.5 ±0.1 






Z(Z ) 


0.24±0.05 






Lx 


2.0±0.1 






X 10 43 erg s _1 








N H xl0 20 cm " 3 


3.3±0.2 




X 2 (d.o.f) 




235 (235) 



NOTES: 

Luminosities quoted are absorbed luminosities in the energy 
range 0.5-5.0 keV 




channel energy (keV) 



Figure 6. The spectrum and residuals to the best-fitting model of source 1 , 
fitted with a MeKaL model. 



Table 6. The best-fitting parameters for the galaxies, fitting a MeKaL 
model plus power law. The absorption was fixed to the Galactic value, and 
the metal abundance of the MeKaL model to 0.3Zq. 



Source 


Component 


Parameter 




A 160 


MEKAL 


kT ( keV ) 
N H (cm -3 ) 


1.23 ±0.4 
4.3 xlO 20 




Power Law 


r 

Lx 

Xl0 41 (erg s" 1 ) 


2.22±0.27 
(1.7±0.3) 




X 2 (d.o.f) 




16.5 (15) 


A 2462 


MEKAL 


kT ( keV ) 
N H (cm" 3 ) 


°- 79 ±o.i 
3.11 xlO 20 




Power Law 


r 

Lx 

XlO 41 erg s - 1 


1.9±0.2 
(5.8±0.2) 




X 2 (d.o.f) 




10.4(13) 



NOTES: 



Luminosities quoted are absorbed luminosities in the energy 
range 0.5-5.0 keV 



2000) (c.f. 2=0.0447 for the cluster), and there is at least one spec- 
troscopically confirmed 'companion' galaxy at a similar redshift 
(z=0.0470; Pinknev et al. 2000). Hardness maps of this region (see 
Section l431 show no evidence of a bow shock in front of the galax- 
ies, and the extended X-ray halo suggests that the galaxies have 
not been significantly stripped by any encounter with the cluster. 
Thus, there appears to be little evidence that this source is interact- 
ing with the cluster in any significant way, and we conclude that it 
is unimportant to the current discussion. 

4.3 Deprojection analysis 

As mentioned in Section Q it is suggested that the external envi- 
ronment is a significant factor in causing jets in WATs to flare sud- 
denly. Therefore, it is important to determine how the properties of 
the cluster vary as a function of position. We constructed tempera- 
ture and density profiles of the ICM across the cluster to see if there 
are any changes in environment with distance. 
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Figure 8. Schematic diagram of an onion skin deprojection. We view the 
cluster from the right hand side; the dotted lines show the line of sight. The 
cluster is modelled as a set of concentric spheres, labelled 1,2. . . n in the 
diagram. The hatched area of the first sphere contributes to some of the 
emission we see from the annulus corresponding to the second sphere. 



For each source, spectra were extracted in variable-width an- 
nuli centred on the radio core, such that each annulus contained 
at least 1000 counts. In Abell 160, we extracted spectra out to 
5.0 arcmin, and in Abell 2462, out to 3.0 arcmin. For Abell 160, 
background regions consisted of scaled blank-sky spectra, taken in 
the same regions as the source spectra. For Abell 2462, we chose 
a circular annulus between 3-4 arcmin, and used that as a back- 
ground for all spectra. Each spectrum was then fitted in XSPEC 
with a MeKaL model, modified by the Galactic absorption. The 
metal abundances and Galactic hydrogen column density values 
were fixed to the values found in Section l4~2l above. The projected 
temperature profiles are shown in Fig. From the profiles, it is 
clear that both clusters have cool cores, which are on similar length 
scales as that of the WAT jets. 

The data were then deprojected, using an onion skin depro- 
jection method. In an onion-skin deprojection, the cluster is mod- 
elled as a series of concentric spheres, as showing in Fig. [8] The 
outermost sphere is first fitted with a MeKaL model, with the 
temperature and normalization free to fit the data. This determines 
the temperature and normalization that is appropriate for the out- 
ermost sphere (labelled 1 in Fig. |8j. The hatched area of sphere 
1 contributes to some of the emission observed from the next an- 
nulus, whose radii correspond to those of the sphere 2 in Fig. [8] 
The MeKaL normalization is volume dependent. The normaliza- 
tion of the first sphere is rescaled to reflect the volume that con- 
tributes to the emission from the second sphere. This first, renor- 
malized MeKaL component is called component 1. The second 
sphere is then fitted with the frozen component 1, and a second 
MeKaL component whose temperature and normalization are left 
free to fit the emission from the second sphere. This second com- 
ponent is renormalized in a similar way to component 1, and the 
frozen component 2 and component 1 are then fitted, together with 
a third MeKaL component to the third sphere. This process contin- 
ues until the (n-l)" 1 sphere has been fitted. The innermost sphere 
is then fitted with the frozen, rescaled, n-l MeKaL components, 
a final MeKaL component, and in our case, a power law compo- 
nent in order to take account of emission from the active nucleus. 
The temperature and normalization of the MeKaL component, and 
the normalization of the power law are left free to fit the data. The 
power law index was fixed to the value obtained in Section l4~2l 

The deprojected profile is presented with the projected tem- 
perature profiles, in Fig.Q There appear to be no major differences 
between the shapes of the two profiles which both agree to within 
la. Both show a steep temperature gradient in the region where the 



jet flares, with the profile then showing only a very shallow temper- 
ature gradient as we move away from the cluster centre. 
The normalization of the MeKaL model is given by 

10" 14 f 

Nmek = : ; -to / n e n p dV (2) 

Atv(D a (1 + z)) 2 J 

where N mB k is the MeKaL normalization for a particular annulus, 
Da is the angular diameter distance of the source, and n e and n v 
are the electron and proton densities respectively. We further as- 
sume that n v = 1.18n e , and that the MeKaL component is fitted 
in a volume V. This implies that 

_ f 4-- [Da (l + z)] 2 N mek \^ 
ne -\ 1.18 x 10---V J ( ) 

The electron density profiles which are obtained in this man- 
ner from the deprojected temperature profiles are shown in Fig. [9] 
From Fig.|5| it can be seen in both profiles that the density decreases 
with distance from the central galaxy. In Abell 160, we see a rather 
steep density gradient over the initial 50 arcsec. This then flattens 
out, and decreases at a slower rate, out to 350 arcsec. The change 
in density gradient does not appear to coincide with the jet flar- 
ing point, and in fact occurs some distance after the point at which 
the jets flare. The second cluster, Abell 2462, has a smooth density 
gradient which does not show any sudden change; the density drops 
rather gradually across the cluster, and the density gradient appears 
smooth at the jet flaring point (at 28 arcsec). 

This discussion, however, does not take into account projec- 
tion effects. We assume here that the radio sources lie on the 
plane of the sky; however, if they were inclined at some angle 
6 to the line of sight, then we would expect the jets to appear 
shorter than they actually are. We can make some estimate of 8 
by considering the properties of the radio jets. The speed of fluid 
flow through the jets is thought to be mildly relativistic, with a 
sp eed o f approximatel y 0.3c (es timated from jet/counter-jet ratios 
by Hardcastle & Sakelliou 2004). Therefore, we should see the ef- 
fects of relativistic Doppler boosting in the jets; the jet travelling 
towards us (the 'jet') will be significantly brighter than the jet trav- 
elling away from us (the 'counter-jet'). Using the radio maps, we 
findjet/counter-jet ratios of 1.1 for 01 10+152 and 2.3 for 2236-176. 
These give inclination angles of 3.6 and 33 deg for 0110+152 and 
2236-176 respectively, which imply that the deprojected jet lengths 
would be 40 arcsec for 0110+152 and 33 arcsec for 2236-176. As 
expected, the projection effects are not great, since the jets and 
counter-jets in both sources are of approximately equal brightness, 
and so projection does not drastically change any of our conclu- 
sions here. 

To obtain the pressure profiles, we treat the plasma as an ideal 
gas and calculate the pressure using P = 2.96 x 10~ w n e kT. Here, 
n e is the electron number density calculated above in cm ~ 3 , kT 
is the temperature of the cluster gas in keV, and we assume that 
n ~ 1.8n e . The pressure profiles derived from the deprojected 
density and temperature profiles are shown in Fig. llQI We note that 
there is a pressure drop across the cluster, as would be expected. 
Also plotted on Fig. ^| are the minimum pressures calculated at 
the bases of the radio lobes of each source. We assume power law 
electron spectra with electron Lorentz factors between 1 and 2 x 10 4 
with the power-law index of the energy distribution equal to -2. We 
also assume that the sources are close to the plane of the sky (as 
confirmed by our calculations above), and cylindrically symmet- 
ric geometry for the plumes. These minimum pressures are signifi- 
cantly lower than the external thermal pressures, by a factor similar 
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Figure 7. The temperature profiles centred on the galaxies - (a) is the profile of Abell 160, and (b) is the profile of Abell 2462. The metalicities and hydrogen 
absorbing column densities are set to the values found for the cluster in Section l4~2l There appear to be definite temperature gradients across the galaxy and 
cluster. The closed points are the projected profiles and the open points are the deprojected profiles. The dashed vertical lines show the jet termination point 
for each radio source. 
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Figure 9. The deprojected electron density profiles for the two clusters - graph (a) is Abell 160 and graph (b) is Abell 2462. 



to those seen in the lobes an d plumes of more typical FRI sources 
lHardcastle & Worrall200d) . 

The temperature profiles in Fig.0 and the surface brightness 
profiles of Fig.[3]are used to define a galaxy-cluster interface. From 
Fig.0 it is seen that galactic emission (at a significantly lower tem- 
perature than the extended emission, as shown by Fig. [7} domi- 
nates the brightness profile out to 40 and 20 arcsec for Abell 160 
and Abell 2462 respectively. Fig. Q shows that the two innermost 
annuli of both clusters have significantly lower temperatures than 
the outer bins. Together, these profiles imply that on a scale of ap- 



proximately 20-40 arcsec, there is an 'interface' between where the 
galaxy dominates, and where the cluster dominates, with a smooth 
transition between the two. 

We do not see any evidence for a shock at the jet flaring point; 
even though the temperature increases rapidly, there is no corre- 
sponding jump in the density gradient at that point. Rather, the 
density continues to decrease steadily, which is counter-indicative 
of a shock. To investigate this in more detail, we also created adap- 
tively binned hardness maps of the central regions of both clusters 
to search for any variation in hardness that would be indicative of 
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Figure 10. The pressure profiles for the two clusters (Abell 160 is on the left, and Abell 2462 is on the right). The vertical lines show the jet flaring position 
for each source. The dotted and dash-dotted lines show the minimum pressure, over the distance in which the pressure was calculated, at the base of the North 
and South lobes of 01 10+152 in the plot for Abell 160, and for the East and West lobes in 2236-176 in the Abell 2462 plot. 



a strong shock that would disrupt the jet. We used a soft energy 
range of 0.5-1.5 keV and our 'hard' energy band was from 1.5- 
3.0 keV . The maps were made usin g an adaptive binning process 
detailed in Sanders & Fabian (2001) to create hardness maps and 
error maps. Overlaying contours of the radio sources (Fig.[2j onto 
the temperature maps shows no significant changes in hardness in 
the regions where the jets flare, compared to regions near to the jet 
flaring point. Therefore, it seems that no strong shocks are seen in 
the data. Some weak shocks may be present, but it is unlikely that 
such a weak shock could disrupt jets such as these, which are likely 
to be moving at speeds much greater than the sound speed in the 
ICM (~0.3c for the jet velocity, compared to ~ 5 x 10 5 ms _1 for 
the sound speed in the cluster). 



5 DISCUSSION 

From the spatial analysis of the X-ray data (Section |4.1> . the jets 
in both sources flare on similar length scales to that of the galaxy- 
cluster transition. It appears from our analysis of the X-ray data that 
most existing models fail to account for the jet flaring in WATs. 
It is known that the environment must be playing a role in deter- 
mining the jet properties, but it is not clear how. Looking at the 
surface brightness and density profiles (Figs.[3]and the density 
appears to vary smoothly across the point where the jets flare in 
both sources. The X-ray gas temperature profiles of both sources 
show a significant increase in temperature across the point where 
the jets flare. We suspect that this is an unresolved steep temper- 
ature gradient, particularly as Chandra observations of the WAT 
3C465 I Hardcastle et al 2005, in prep.) have sufficient signal to 
noise to detect a temperature gradient in the equivalent region of 
that source. In all these sources, and in archival observations of 
Hydra A (which we have reduced in a similar way to the data pre- 
sented here) the striking feature is that the transition between the 



inner, well-collimated jets and the diffuse plumes occurs on a very 
similar size scale to the transition between a cool, dense 'core' at 
the centre of the cluster and the hotter, more diffuse cluster emis- 
sion. This leads us to believe that the cool core may be closely 
related to the location of the jet-plume transition. We discuss here 
models put forward to explain the jet flaring process. 

5.1 Models unable to explain the jet flaring 

• The lLoken et alj 1 19951) 'cross-wind' model, as detailed in 
Section Q is discussed in lHardcastld Il99^) in relation to 3C130. 
The arguments put forward against the cross-wind model in that 
paper also apply to 01 10+152 and 2236-176; bent jets are not seen, 
and neither do the plumes bend immediately after the jet disruption 
point, as th e cross-win d model would imply. 

• lLoken et al] ll995) also put forward a second model to explain 
jet disruption - they model the jet as disrupting when it propagates 
across a shock in the ICM. We would expect to see some evidence 
of shocks in our data, if they are responsible for causing the WAT 
jets to flare. However, there is no evidence for a shock in our Chan- 
dra data (see Section l4~3l . and so this model also fails to explain 
the j et disruption. 

• iHiggins et alj J 19991) model a jet passing through a medium 
containing cool dense clouds, which are denser than the ambient 
medium by a factor of 50, but are in pressure balance with the ICM. 
They model relativistic jets, in pressure balance with the ICM, hit- 
ting these dense clouds. The model predicts that moderately fast 
jets should disrupt and flare suddenly on coming into contact with a 
cloud, and that the radio source should resemble a WAT. We exam- 
ined our X-ray images for evidence of cool clouds spatially coinci- 
dent with the jet flaring points. In both sources, we find no evidence 
for any such clouds. Furthermore, finding clouds equidistant from 
the host galaxy, which disrupt the jets seems rather implausible - 
clusters are large when compared to clouds and jets, and having 



10 N. N. Jetha et al. 



We consider the radial component of the Navier-Stokes equa- 




tion: 



Figure 11. A schematic diagram of what may be happening in a WAT 



clouds placed equidistant from the host galaxy in all WAT clusters 
seems too far-fetched a coin cidence. 

• iHooda & Wiital Il996) model a jet passing through a power 
law ISM, and subsequently a constant density ICM. The jet be- 
comes unstable to the growth of Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities on the surface of the jet and the jet co- 
coon upon crossing the ISM/ICM boundary. Their simulations 
find that low Mach number jets disrupt on crossing the bound- 
ary, and furthermore, disrupt faster than their high Mach number 
counterparts. Whilst their simulations produce flaring jets, they as- 
sume abrupt gradients in the density and temerapture of the ex- 
ternal medium that are not consistent with our observations. Also, 
given the best estimates of W AT jet speeds, (approximately 0.3c; 
lHardcastle & Sakellioull2004) . it is likely that WAT jets have high 
Mach numbers with respect to the ICM. To fully understand the 
processes involved in forming WATs, realistic simulations of radio 
galaxies in clusters are needed; such simulations need to include 
models for the ICM that match the data (i.e. continuous density 
and steep temperature gradients separating the cool core from the 
rest of the cluster) as well as simulating radio jets with higher Mach 
numbers. 



5.2 An Alternative Model? 

Up until recently, it has been assumed that the hot X-ray emitting 
ICM is responsible for the jet termination. However, as mentioned 
previously, in almost all cases studied, the jet terminates inside the 
radio plume. Therefore, we may have been misled in thinking that 
the jet termination is determined by the ICM properties. It might be 
more appropriate to assume that the the conditions inside the radio 
plume cause teh jet termination, and the external ICM determines 
properties such as the shape and length of the radio plumes. 

We look to simulations of radio galaxies in clusters for ver- 
ification of this. Om ma et alj 12004) study how radio galaxies in 
systems with cool cores heat the ICM, and how the radio galaxies 
themselves develop over time. Their simulations of a cluster con- 
taining a central radio galaxy initially show a source akin to a small 
FRII, but as time progresses, this radio galaxy gradually starts to 
look like a WAT. The jet material then flows back from where the 
jet terminates towards the cluster centre (their fig. 4). How far to- 
wards the cluster centre this flow extends to must be determined by 
where the pressure gradient in the cluster halts the jet material. 

We consider the situation in Fig. II II material leaves the end 
of the jet, at a distance r, e t> with a velocity v max , and under the 
influence of gravity flows back up the pressure gradient to the point 
Tmin where it is stopped by the pressure gradient. This should give 
us a theoretical estimate as to where the base of the plume should 
be located. 



Pplumc^ ~ 



dv 



dPplun 



GM(< r) 



Ppiumc 



(4) 



dr dr r 

here, p p i U mc is the density of the radio plume, P p i umc is the internal 
plume pressure, M(< r) is the mass of the cluster within a radius 
r from the centre of the cluster. 

We then assume that the cluster is in hydrostatic equilib- 
rium, and that the lobe and ICM are in pressure equilibrium, so 
dPpiume/dr = dP icm /dr, and that p p i ume < picm, and obtain: 

d« 2 1 dP ic 

~ -2 

dr ppiumc dr 

Integrating both sides from where the jet terminates (at rj ct , v 
i>max) to where v = at r m t n , we get: 



(5) 



1 dP ic 



Ppiumc 



dr 



-dr 



(6) 



In order to apply this model to the real data, we can 
fit the observed densities and pressures, from Figs. [9] and 1101 
with power laws to simply model the state of the ICM. Using 
Abell 160 as an example, we obtain p cm = Po (r/r ro f )~ - 25 
and picm = po (r/r ro f ) -0,84 , where Po = 5.35 x 10 _13 Pa, and 
po=4130amu m -3 . 

The distance from the radio core to the base of the plume is 
~35 kpc in this source. This can be reproduced from the fits to 
the data and the above model by taking suitable values for w max 
and Ppiumc whic h in this case are n ma x = lOOOkms" 1 (from 
Omma et al. 2004) and p p i um c = 0.001pi C m- The value for p p i um e 
is, however, several orders of magnitude larger than what we would 
expect for a fully relativistic plasma, but given that the lobe could 
be entraining colder gas from its environment, p p i um c = 0.001pi C m 
is not entirely unrealistic. 

It appears from the preceeding discussion that we could turn 
a weak FRII into a WAT by placing it in a cluster with a steep 
pressure gradient near the core, to surpress the back-flow. This 
implies that there should be some sort of relationship between 
how far upstream the back-flow travels and the temperature of 
the cluster. To this en d , we use the WAT sample presented in 
lHardcastle & SakellioJ 120041) . to determine if there exists any 
correlation between rj e t — r m in and the cluster temperature. We 
present our findings in Fig.^| which shows an anti-correlation be- 
tween rjct — r m in and the cluster temperature. This suggests that 
the higher pressures in the centres of hotter clusters may suppress 
the back-flow. This may result in conditions in the radio plumes 
of hotter clusters being altered so that jets in those clusters prop- 
agate through smaller distances in the radio plume before being 
disrupted. 

From our results, it appears that a cool core, joined to the 
rest of the cluster by a region with a steep temperature gradient, 
is an essential feature of any WAT host cluster. It is not clear how 
such a cool component might cause a WAT to form from a weak 
FRII. One possibility is that the link arises from the higher pressure 
within cool cores. From hydrostatic equilibrium it can be seen that 
for a given value of M(< r)/r 2 , the pressure gradient is steeper 
when picm is larger. Hence, in dense cool cores, the steeper pressure 
gradient could inhibit back-flow more effectively, giving a plume 
whose bas e is closer to the je t termi nation point. 

Hard castle & Sakelliod l2004h also find a significant anti- 
correlation between rj c t and the cluster temperature - jets in hotter 
clusters are shorter than those in cooler clusters. This suggests that 
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the Hardcastle & Sakelliou 1 2004) sample of WATs. There appears to be 
an anti-correlation between the cluster temperature and rj ct — r m ; n 

the cluster is a factor in determining the location of the jet termina- 
tion. The trend for shorter jets in hotter systems may be related to 
their higher gas pressures, but simulations of weak FRIIs in clusters 
that contain realistic cool cores are needed to investigate this. 



6 CONCLUDING REMARKS 

We find that new Chandra data shows a relationship between the 
radio morphology of WATs, specifically the jet length, and the tem- 
perature distribution of the cluster medium. In particular, it appears 
that the location of the base of the plume is correlated with the steep 
temperature gradient that separates a cool core from the rest of the 
ICM. However, we are still unclear on the mechanism behind jet 
flaring and the locating of the base of the plume. Firstly, we need to 
establish whether a cool core on scales similar to the length of the 
radio jets is the exception rather than the rule; secondly we need to 
understand why the steep temperature or pressure gradient causes 
jets to flare. The first point requires deeper X-ray observations of a 
large, statistically significant sample of WAT host clusters to com- 
pare with both clusters that host radio quiet galaxies, and clusters 
that host 'normal' FRI galaxies. For the second point, the properties 
of WAT jets and of their cluster environments need to be fed into 
realistic cluster models in order to establish whether the steep tem- 
perature gradient is in itself responsible for the jet flaring, or rather 
is simply a tracer of some underlying phenomenon that causes the 
jets to make the sudden transition to a diffuse plume. 
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